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Nonlinear integrable discrete systems have recently attracted much at- 
tention. Among them, the most fundamental ones are the discrete analogues 
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of Painleve equations. Some properties such as the existence of Lax pair, 
Backlund transformation and singularity confinement property, have been 
studied for these discrete Painleve equations |I|. However as for the solu- 
tions, only in a few cases particular solutions have been derived and the 
structure of the solution space is still almost unknown. We have studied the 
solution in determinant form for the discrete Painleve II equation based on 
the bilinear formalism^]. We here investigate the bilinear structure for the 
discrete Painleve I (dPI) equation on semi-infinite lattice, and show that its 
solutions are given in terms of the Casorati determinant. 

The dPI equation on semi-infinite lattice, 

bN 
u N+1 + u N + u N -i = a + - — , for N > 1, u = 0, (1) 

2uj\[ 

where a and b are arbitrary constants, appears in the theory of two-dimensional 
quantum gravity 0. The partition function of the matrix model, 

Z = fdMexp(-trV(M)), (2) 

V(M) = ciM 2 + c 2 M\ (3) 

where the integral is taken over all m X m Hermite matrix M with the Haar 
measure, is given in terms of the solution of dPI eq.(l). Actually we have 



Z = m\w m u 1 m ~ 1 u 2 m - 2 ---u m - U w= dxexp(-V(x)), 

J — oo 

where u^ is the solution of (1) with the boundary condition, 

1 r°° 
U\ — — I dxE 2 exp(— V(x)), 

W J-oo 

and a, b and ci, c 2 are related as 

2c 2 ' 2c 2 ' 

We expect that investigation of solutions of discrete Painleve equations is 
useful to understand the structure of physical systems such as the above 
matrix model. 



EXACT SOLUTIONS OF DISCRETE PAINLEVE I EQUATION 



In the context of the bilinear formalism, the dPI eq.(l) is naturally ex- 
tended to the discrete coupled system, 



v n N + u n N + v n N _ l = a+^- 1 forJV>l, 

u n N+l + v n N + u n N = a+ b{ - n \P + C , foriV>0, 



(4) 



V N 



with the boundary condition, 



< = 0. 



In (4), n is the auxiliary suffix introduced to have a natural expression for 
the solutions in bilinear formalism. The coupled system (4) reduces to the 
dPI eq.(l) by taking c = 6/2, n = and 

U2N = u%, U 2N +1 = v%, for N >0. 

In the following, we investigate the solution of this extended dPI equations 
using the bilinear method. The bilinear method is originally developed to 
solve the soliton equations^]. Recently it has been shown that the method 
is also applicable to the discrete Painleve II equation^]. 

First we introduce the r function by the dependent variable transforma- 
tion, 

n _n+l _n+l _n 

,n __ 'N+l' N-1 „,n _.. T N+1'N 



Uat — .it _ , V 



v -n+1 n ' U N n n +l ' 

r JV 'N 'N+l'N 



Substituting these into (4) and decoupling the resultant equations, we arrive 
at the bilinear form for the extended dPI equations, 

_n— 1 n+l n-1 n+l n n ( r\ 

T N+1 T N-1 — T N T N T N T Ni \°) 

t n +iT n _ x = (Kn + N-l)+ c)r n N r n N - (b(n - 1) + c^+V^S (6) 

T N+1 T N ~ T N+2 T N-1 = aT N+l T N + {K n ~ *) + c ) T N+l T N • {') 

Equations (4) actually follows from (5)- (7). For example we derive the first 
equation in (4). Dividing (7) by t n+1 t n , we get 

n+l n-1 n-1 n+l n-1 n+l 

Tn+iTn - Tjv+2Tjv - 1 = a + (bin - 1) + c ) Tn+iTn , 

r N+l T N r N+l r N T N+l r N 
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which is rewritten by using (5) as 



v"% + u n N = a + (b(n - I) + c) 



ra-l n+l 

T N+1'N 

T N+1'N 



(8) 



Eliminating the term rjy t n from (5) and (6), we have 



T N+1 T N-1 



hNr n T n 



{b{n - 1) + c)t 



n— 1 _n+l 
TV+l'TV-l- 



Multiplying this by r" +1 /(r^ +1 r^r™+ 1 1 ), we get 

hN T n-1 n+1 

«7V 



r iV+l r iV 



(9) 



Adding (8) and (9), we obtain the first equation in (4). The second equation 
in (4) is similarly derived from (5)-(7). Finally the boundary condition Uq = 
is satisfied by taking r" x = 0. 

Now we solve the bilinear difference equations (5)- (7). The solution is 
given as 



'N 



A 



A n 


A n +i 




A n+ N-l 


A n +\ 


A n +2 




A n +N 


-n+N- 


L A n+N 




A n +2N-2 




T n — 1 

'o ~~ -"-> 


t\ = 0, 



N>h 



(10) 



where A n is a discrete Airy function satisfying 

A n+2 = aA n+1 + (bn + c)A n . 



(11) 



It is proved by using the standard technique || that the above r^ actually 
satisfies (5)-(7). As an example, we give the proof for (6). In the Casorati 
determinant (10), adding (N — l)-th column multiplied by —a to iV-th col- 
umn, adding (N — 2)-th column multiplied by —b(N — 2) to iV-th column 
and using (11), we get 



t 



N 



A 



A n 


A n +l 


A n+ N-2 


(b(n - 1) + c)A n+N _ 3 


A n +\ 


A n +2 


An+N-l 


(bn + c)A n+N ^ 2 


n+N-l 


A n +N 


' A n+ 2N-3 


(b(n + N-2)+ c)A n+2T , 
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Furthermore, adding (N — j — l)-th column multiplied by —a to (N — j)-th 
column and adding (N — j — 2)-th column multiplied by —b(N — j — 2) to 
(N — j)-th column, for j = 1, • • ■ , N — 2 successively, we have 



t. 



A' 



(b(n - 1) + c)(bn + c) ■ ■ ■ (b(n + N - 2) + c) 



i A 

b{n-l)+c n 
ta+c n+1 



A n -! 

A„ 



A n +N-3 
A n +N-2 



b ( n+N -2)+c A ™+N-l Ai+N-2 



A 



n+2N-A 



Applying the Jacobi formula to the above determinant, we obtain the bilinear 
equation, 



T 



N 



(b(n - 1) + c)(bn + c) • • • (b(n + N - 2) + c) 



<N-2 



T, 



N-l 



{b{n - 1) + c)(bn + c) ■ ■ ■ (b(n + N - 3) + c) 



'7V-1 



T 



n+1 



AT-1 



(bn + c)(b(n + 1) + c) • • • (b(n + N - 2) + c) 



r 



n-l 
JV-D 



which is essentially the same as (6). Hence we have shown that r^ given by 
(10) satisfies the bilinear difference equation (6). The other equations are 
also proved in a similar way. 

We note that the discrete Airy function defined by (11) is expressed as 

A n -{ 0) iln-1+c/b I /_ b y/2 j ' 

where H u (x) is the Hermite- Weber function satisfying the contiguity relation, 
H u+ i(x) - xH v (x) + vH v _i(x) = 0. 



We have investigated the Casorati determinant solution of the extended 
dPI equations. As was mentioned, the solution for dPI equation gives the 
partition function of a matrix model. Hence a natural question is what kind 
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of matrix model corresponds to the extended dPI eqs. (4). We find that it 
is the matrix model with the potential, 

V(M) = ciM 2 + c 2 M 4 + c logM 2 . (12) 

Actually the partition function defined in (2) with the potential (12) is given 
as 

Z = m\w m (v° ) m ~\u° l ) m - 2 {v i i) m -\u Q 2 ) m ~ 4 

where u° N and v% are the solution of (4) with the boundary condition, 

1 f°° 
Vn = — dxx 2 exp(— V(x)), 

W J-oo 



and the parameters, 



c\ _ 1 _ 1 -2c 

"2c 2 ' ~2c 2 ' ( '~ 4c 2 



Finally we comment on the continuous limit of dPI. The dPI equation 
itself reduces to the Painleve I equation on some condition for the parameters. 
However, the relation between the solution given in this article and that 
for the Painleve I equation is not clear yet, because the r function (10) is 
considered only on the semi-infinite lattice. 
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